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Abstract 

We suppose that there is new physics (NP) in b — > s transitions, and examine its 
effect on the angular distribution of B® — > ViV 2 (q = d,s), where V\ t i are vector 
mesons. We find that, in the presence of such NP, the formulae relating the param- 
eters of the untagged, time-integrated angular distribution to certain observables 
(polarization fractions, CP-violating triple-product asymmetries, CP-conserving in- 
terference term) must be modified from their standard- mo del forms. This modifica- 
tion is due in part to a nonzero B®-B® width difference, which is significant only for 
B° s decays. We re-analyze the B® — > <fxj) data to see the effect of these modifications. 
As AT S /2T S ~ 10%, there are 0(10%) changes in the derived observables. These 
are not large, but may be important given that one is looking for signals of NP. In 
addition, if the NP contributes to the b — > s decay, the measurement of the untagged 
time-dependent angular distribution provides enough information to extract all the 
NP parameters. 
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1 Introduction 



As recently as a year ago, there were several hints of physics beyond the stan- 
dard model (SM) in b — > s transitions. For example, the CDF pQ and D0 [2] 
Collaborations measured the CP asymmetry in B° s — > J/ip(j), and found a hint for 
indirect (mixing-induced) CP violation. This is counter to the expectation of the 
SM, which predicts this CP violation to be ~ 0. In general, this result was inter- 
preted as evidence for a nonzero value of the weak phase of B Q S -B® mixing (2/3 s ), 
and the contributions of various new-physics (NP) models to the B s mixing phase 
were explored [3], SI El El d [HI [9] . It was also pointed out that NP in the decay 
b sec could also play an important role [10]. In addition, the SM predicts that 
the measured indirect CP asymmetry in b — > sss penguin decays should gener- 
ally be equal to that found in charmonium decays such as B® — > J/ipKs- How- 
ever, it was found that these two quantities were not identical for several decays 
[TT] . As a third example, the CDF Collaboration reported the measurement of 
B(B° S = (l-Stol) x 10 " 8 H3- This is lar § er than the SM prediction for 

this branching ratio, which is B(B° S = (3.35 ± 0.32) x 10~ 9 [33]. There 

were a number of other effects - in all cases, the disagreement with the SM was not 
large, < 2a. Still, it was intriguing that all appear in b — > s transitions. 

However, with recent measurements, these effects have largely disappeared, or 
at least been reduced. First, LHCb has measured the indirect CP asymmetry in 
B® — > J/ip4>, and finds 2(3 s ~ 0, in agreement with the SM [13] . Specifically, they 
measure 2(3 S = (—0.06 ±5.77 (stat) ± 1.54 (syst))°. Still, the errors are large enough 
that NP cannot be excluded. Second, with the latest indirect CP asymmetry data, 
the Heavy Flavor Averaging Group (HFAG) [T5] finds that the B^-B^ mixing phase 
sin 2/3 is measured to be (i) 0.68 ± 0.02 in charmonium decays, and (ii) 0.64 ± 0.04 
in b — > sss penguin decays. These numbers are quite similar, so that no real 
discrepancy can be claimed. On the other hand, several of the b — > sss decays 
have additional contributions with a different weak phase, and so HFAG urges that 
the "naive average" in (ii) be used with extreme caution. Third, the recent LHCb 
update does not confirm the CDF B° s — > fi + n~ result [16] . They improve the present 
upper bound to B(B Q S < 1.3 x 10~ 8 (90% C.L.), in agreement with the SM. 

Most of the other effects have similarly gone away, or are simply not large enough 
to be compelling. 

On the other hand, there is one discrepancy with the SM which has not dis- 
appeared. The D0 Collaboration has reported an anomalously large CP-violating 
like-sign dimuon charge asymmetry in the B system. In Ref. [17], the asymmetry 
was found to be 

A b sl = -(9.57 ±2.51 ± 1.46) x 10~ 3 , (1) 
which is a 3.2a deviation from the SM prediction, A b f M = (— 2.3±g x 10" 4 [H]. 
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In fact, the updated measurement [19] exhibits an even larger discrepancy: 

A b sl = -(7.87 ± 1.72 ±0.93) x 1(T 3 , (2) 

a 3.9cr deviation. Now, it has been shown that this anomaly is due mainly to B s 
decays, i.e. aft-ls transition. So this is a solid indication of b — > s NP. 

Thus, at present the status of b — > s NP is uncertain. It seems unlikely that the 
effect of such NP can be very large, but a smaller effect is still possible. In this paper 
we make the assumption that NP is present in b — > s transitions. However, in addi- 
tion to taking into account its effect on B s mixing, which is what is conventionally 
done, we also consider its effect on b — > s decays. The main aim is to examine the 
effect of b — > s NP on the angular distribution of B® — > V1V2 (q = d,s), where 
are vector mesons. In particular, we consider final states which are self-conjugate, 
so that both B® and B® can decay to V1V2, generating indirect (mixing- induced) 
CP-violating effects. 

There are three classes of B® decays which can be affected by b — > s NP: 

1. B® decays with b — > s, 

2. B® decays with b — > s, 

3. B° s decays with b — >■ d. 

Our analysis is completely general and can be applied to any of these classes. How- 
ever, we also focus specifically on B® — > cfxf). There are two reasons. First, this is a 
pure b — > s penguin decay, and so there can well be NP contributions to any of the 
loop-level penguin decay amplitude^]. Second, the untagged angular distribution 
of the decay has already been measured by the CDF [21] and LHCb [22] Collab- 
orations, and so their results can be (re) interpreted in the context of b — > s NP 
contributions. 

The result of this analysis - and this is the main point of the paper - is as 
follows. The parameters of the untagged, time- integrated angular distribution can be 
measured experimentally. Certain observables can be derived from these parameters. 
However, in the presence of NP, the formulae which relate the observables and 
parameters are modified compared to their SM forms. There are six terms (i = 1-6) 
in the angular distribution, and we correspondingly find six observables for which 
the relation between the experimental data and theoretical parameters must be 
modified. For i = 1-3 they are the polarization fractions, for i = 4,6 they are the 
CP-violating triple-product asymmetries, and i = 5 corresponds to a CP-conserving 
observable. The modifications for the polarization fractions are particularly striking. 
Here there are corrections to the SM formulae that are proportional to the width 
difference in the B®-B® system. Now, the width difference Ar is sizeable only for B° s 

— > (jxf) and — > J/ip(f> were examined in Ref. [20], but only NP in B®-B® mixing was 
considered, not NP in the decay. 
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decays. Thus, the formulae modifications due to NP are important only for class-(l) 
and (3) decays, which include B® — > cfxp. Ar s /2r s ~ 10%, so that the modifications 
lead to 0(10%) changes in the derived observables. These are not large, but may 
be important given that one is looking for signals of NP. 

Another result is that, if the untagged, time-dependent angular distribution can 
be measured, 12 observables can be obtained. If the NP contributes to the b — > s de- 
cay, there are fewer than 12 unknown NP parameters. Thus, all of these parameters 
can be extracted from the angular distribution. This may allow the identification of 
the NP. 

We begin in Sec. 2 by presenting the most general B® s — > V1V2 angular distri- 
bution, allowing for NP in the mixing and/or the decay. We consider the angular 
distribution for several different scenarios: at t — (Sec. 2.2.1), time-dependent 
(Sec. 2.2.2), untagged time-dependent (Sec. 2.3), untagged time-integrated (Sec. 2.4). 
In Sec. 3 we examine the untagged time- dependent and time-integrated distributions 
for 5° 00 within the SM. The study of B° s -> 00 is extended to the SM + NP 
in Sec. 4. We discuss observables such as the polarization fractions, CP-violating 
triple-product asymmetries and the CP-conserving interference term, and note the 
changes in the formulae used for their extraction necessitated by the inclusion of 
b — > s NP. We also show that all the unknown NP parameters in the b — > s decay 
can be determined from the measurement of the untagged, time- dependent angular 
distribution. In Sec. 5 we present a numerical reanalysis of the B° s — > 00 data al- 
lowing for the possibility of b — > s NP contributions in the decay. We conclude in 
Sec. 6. 

2 B — ► V1V2 Angular Distribution 

2.1 Generalities 

The most general Lorentz-covariant amplitude for the decay B(p) — > Vi(fci,£i) + 
V 2 (k 2 , e 2 ) is given by [231 [21] 

M = ael-s 2 + \{p ■ el){p ■ e* 2 ) + i-^e^p^e*^ , (3) 

where q = ki — k 2 . The quantities a, b and c are complex and contain in general 
both CP-conserving strong phases and CP- violating weak phases. In B — > V1V2 
decays, the final state can have total spin 0, 1 or 2, which correspond to the V\ 
and V2 having relative orbital angular momentum I — (s wave), I = 1 (p wave), 
or I — 2 (d wave), respectively. The a and b terms correspond to combinations 
of the parity-even s- and (i-wave amplitudes, while the c term corresponds to the 
parity-odd p-wave amplitude. 

In order to obtain the angular distribution for B — > V1V2, one uses the linear 
polarization basis. Here, one decomposes the decay amplitude into components in 
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which the polarizations of the final-state vector mesons are either longitudinal (Aq), 
or transverse to their directions of motion and parallel (A\\) or perpendicular (A±) 
to one another. The transversity amplitudes Ah (h = 0, ||, _L) are related to a, b and 
c of Eq. (j3J) via 



A\\ = V2a, A = -ax - ^^b(x 2 - 1) , A ± = 2^2 ^i^cV^T , (4) 



where x = k\ ■ k 2 / '(mim 2 ) (m\ and m 2 are the masses of V\ and V 2 , respectively.). 

The amplitude for B(p) — > Vi(A;i,£i) + ^2(^2, £2) can be obtained by operating 
on Eq. fl3]) with CP. This yields 

M = a el ■ e* + -A-(p ■ e\){p ■ e* 2 ) - i^-e^^e*^ , (5) 

in which a, b and c are equal to a, b and c, respectively, except that the weak phases 
are of opposite sign. The above equation can be obtained from Eq. ()3]) by changing 
a — > a, b — > b and c — > —c. Similarly, one defines Aq, A\\, and A±, which are equal 
to Aq, A\\, and A±, respectively, but with weak phases of opposite sign. 

2-2 Bl^VM 

As mentioned in the introduction, we are interested in the decays B® — > V1V2 (q = 
d,s), in which both B® and B® can decay to V1V2. Due to B®-B® mixing, the 
amplitude is time dependent. Assuming that Vi )2 both decay into two pseudoscalars, 
i.e. V\ —> P\P[, V2 —> P2P21 t ne angular distribution is given in terms of the vector 
uj = (cos 0i, cos 2 ,$) [211 [26]: 

<M = £± m)m . (6) 



i=l 



Here, 9\ (0 2 ) is the angle between the directions of motion of the Pi (P 2 ) i n the V\ 
(V 2 ) rest frame and the V\ (V2) in the B rest frame, and $ is the angle between 
the normals to the planes defined by P\P[ and P 2 P 2 in the B rest frame. The 
angular-dependent terms are given by 

/1 (3) = 4 cos 2 0i cos 2 2 , f 2 (u) = 2 sin 2 Q x sin 2 2 cos 2 $ , 

f 3 (uj) = 2 sin 2 0i sin 2 2 sin 2 $ , f 4 (uj) = -2 sin 2 Q x sin 2 2 sin 2$ , 

f s (u) = v 7 ^ sin 20i sin 20 2 cos $ , f 6 (uj) = -V2sin20i sin20 2 sin$ . (7) 
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2.2.1 t = 



At t — 0, the Ki are 

K 1 = \A \ 2 , K 2 = |A||| 2 , K 3 = \A ± \\ 
K, = lm(A ± A0 , K 5 = Re(A l{ A*) , K 6 = lm(A ± A*) . (8) 

The angular distribution for the CP-conjugate decay B® — > ViV 2 is the same as that 
given above with the replacements Ki — > Ki, and Ah — > Af L . 

The quantities K4 and K e are particularly interesting. They are related to the 
tfivpaV^I 1 J £\£*2 term of Eq. Q, which is proportional to q- (ei x e 2 ) in the rest frame 
of the B. This is the triple product (TP). The TP is odd under both parity and 
time reversal, and thus constitutes a potential signal of CP violation. However, here 
one has to be a bit careful. As noted above, the Ah possess both weak (CP-odd) 
and strong (CP-even) phases. Thus, K4 and/or K Q can be nonzero even if the weak 
phases vanish. In order to obtain a true signal of CP violation, one has to compare 
the B and B decays. Now, is the same as K^, except that (i) the weak phases 
change sign, and (ii) there is an overall relative minus sign due to the presence of 
Aj_/A±, and similarly for K 6 and K Q . This implies that the true (CP-violating) TP's 
are given by the untagged observables K4 + K4 and K$ + K$. There are also fake 
(CP-conserving) TP's, due only to strong phases of the the AhS, given by K4 — K4 
and Kq — Kq. 



2.2.2 Time dependence 

In order to calculate the Ki(t), one proceeds as follows. Due to B®-B® mixing, the 
time evolution of the states |£>°(t)) and |-B°(t)) can be described by the relations 

\B° q (t)) = g + (t) \B° q ) + q -g4t) \B° q ) , 

\B° q (t)) = - q 9-(t) \B° q )+g + (t) \B° q ) , (9) 

where q/p = e~ l ^ q . Here, (f> q is the phase of B q -B q mixing. In the SM, we have 
tfid = 2(3 = (42.8 ± 1.6)° from charmonium decays [15]. Also, assuming no NP in the 
decay, the LHCb Collaboration measures <fi s = (—0.06 ±5.77 (stat) ±1.54 (syst))° in 
B° s Hlj. Although this agrees with the SM prediction of (j) s ~ 0, the errors 

are still large enough that NP in the decay and/or mixing cannot be excluded. 
In the above, we have 




-{iM L +T L /2)t 1 -{iM H +T H /2)t 
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g_(t) = i( e -(^+ r ^/^ _ e -{iM H+ V H /2)t^ ^ (1Q) 

where L and H indicate the light and heavy states, respectively. The average mass 
and width, and the mass and width differences of the B- meson eigenstates are defined 
by 

_ M H + M L _ T L + r H 

m ~ 2 ' 2 (11) 

Am = M H - M L , Ar = T L - T H . 

Am is positive by definition. For mesons, I\ ~ Th, so that AI^ = 0. How- 
ever, for Bg mesons, AT S is reasonably large: |Ar s | = 0.123 ± 0.029 (stat) ± 
0.011 (syst) ps _1 [28]. In our convention the SM prediction for Ar s is positive, 
and it has been recently confirmed experimentally that Ar s > [29J. 

The time dependence of the transversity amplitudes is due to B®-B® mixing. 
For the decay to a final state / we have 

A h (t) = (f\B° q (t)) h = [g + (t)A h + r) h q/pg-(t) A h j , 

A h (t) = (f\B° q (t)) h = [p/q g_{t)A h + V h g+(t) A h ] , (12) 

where A h = (f\B°) h , A h = (f\B°) h , and 770,11 = 1, V± = -1- In calculating the 
Ki(t), the following relations are useful: 

\g±(t)\ 2 = i e - r *(cosh(Ar/2)t ± cos Ami) , 

g%(t)g-{t) = l -e~ vt ( - sinh (AT/2)t + % sin Ami) . (13) 
The expressions for the time-dependent functions K{(t) are given by 
K x {t) = \A (t)\ 2 = l -e~ vt [{\A Q \ 2 + \A Q \ 2 ) cosh {AV/2)t 

+ (\A Q \ 2 - \ Aq\ 2 ) cos Ami 

— 2He(A* Q Ao) (cos0 q sinh (Ar/2)t — sin</> 9 sin Amt) 

- 2Im(AoA ) (cos <p q sin Amt + sin <p q sinh (AT/2)t)] , 

K 2 {t) = |A||(t)| 2 = X -e~ vt [(|A||| 2 + |A||| 2 ) cosh(Ar/2)t 

+ (|A||| 2 - |A||| 2 ) cos Amt 
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- 2Re(A* ] A\\ ) (cos <j> q sinh (Ar/2)t - sin <p q sin Amt) 

- 2Im(A|A||) (cos <p q sin Amt + sin <p q sinh (Ar/2)t)] , 

A' 3 (t) = |A ± (t)| 2 = V ri [(|AJ 2 + |A ± | 2 ) cosh(Ar/2)t 

+ (|AJ 2 - \A ± \ 2 ) cos Amt 
+ 2Re(y4^_A_L) (cos (f) q sinh ( Ar /2)t - sin (j) q sin Amt) 
+ 2\m(A* L A A _) (cos g sin Amt + sin cj) q sinh ( Ar /2)t)] , 

K A (t) = lm(A ± (t)Al(t)) = X -e~ vt [(lm(A x A|) - Im(A ± A|)) cosh(Ar/2)t 

+ (lm(Aj_A|) + lm(A ± AV)) cos Amt 
+ (lm(A ± v4|) -Im(A ± A|)) (- sinh (Ar/2)t cos 4> q + sin Amt sin <f> q ) 
+ (Re(A±A* ]j ) + Re(A ± Al)) (- sinh (AT/2)t sin - sin Amt cos 0,)] , 

tf 5 (f) = Re(A||(t)AS(t)) = V R [(Re(^||A*) + Re(A,|A*)) cosh (Ar/2)t 

+ (Re(A||A*) - Re(A||A*)) cos Amt 
+ (Re(A|| Aq) + Re(A|| Aq)) (- sinh (Ar/2)t cos 9 + sin Amt sin <f> q ) 
+ (lm(A|| Aq) - lm(A\\AD) (sinh (Ar/2)t sin 9 + sin Amt cos <f> q )] , 

K 6 (t) = lm(A ± (t)A*(t)) = l -e- Tt [(lm(A x A* ) - lm(A ± A*)) cosh (Ar/2)t 

+ (kn{A ± A* ) + lm(A ± A* )) cos Amt (14) 

+ (lm(A ± A* ) - Im( A ± Aq)) (- sinh ( Ar /2)t cos 4> q + sin Amt sin <p q ) 

+ (Re(A ± A*) + Re{A ± A* )) (- sinh (Ar/2)t sin <p q - sin Amt cos </>„)] . 

The expressions for the time-dependent Aj(t)'s can be obtained from the -fTj(t)'s 
by changing the sign of the weak phases in both the decay (Ah i]hAh) and the 
mixing (<p q -)• 
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2.3 Untagged decays 

In the previous subsections, we presented the angular distribution for the case in 
which the initial decay meson is tagged, so that one can distinguish the B® and B® 
decays. In practice, however, tagging is difficult. Thus, as a first step, experiments 
will examine the untagged decay, and this is considered here. 
The untagged time-dependent angular distribution is given by 



d\T B " + T B o 
dtdu 



9 
32^ 



to 



i=l 



where the untagged observables can be found using Eq. ( 1141) : 
Kx(t) + K^t) = e~ Tt (\A Q \ 2 + |A | 2 ) cosh (Ar/2)t 

- 2(Re(A*A ) cos<p q + lm(A* A ) sin^) sinh (AT/2)t 
K 2 (t) + K 2 (t) = e" r *[(j,4||| 2 + |A||| 2 ) cosh(Ar/2)t 

- 2^Re(AJA||)cos0 3 + Im(A|A||)sin^ sinh(Ar/2)t 



K 3 (t) + K 3 (t) = e~ rt ( |AJ 2 + \A ± \ 2 ) cosh (Ar/2)t 



(15) 



2 Re( A^ A ± ) cos <p q + Im( A^ A ± ) sin <p q sinh ( AT/2)t 



iiT 4 (t) + Ki{t) = e~ rt ( Im(A ± A|) - Im(A ± A|) ) cosh (Ar/2)t 



(Im(A ± A|) - Im(A ± A|)) cos, 



+ (Re(A ± A|) + Re(A ± AjP) sin0J sinh (Ar/2)t 



K 5 (t) + lf 5 (t) = e~ rt (Re(A||A*) + Re(A {l A* )) cosh (AT/2)t 



- ((Re(A||A*) + Re(A||A*))co Sl 



- (Im(A||A*) - Im(A||A*)] sin g sinh (Ar/2)t 
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K 6 (t) + K 6 (t) 




(lm(A ± A* ) - lm(A ± A* )) cos0 g 



+ (Re{A ± A* Q ) + Re(A ± A* Q )) sin^J sinh (Ar/2)t 



(16) 



Note that the CP properties of all the terms are respected. For example, the 
Ki(t) + Ki(t) (i = 1,2,3,5) are supposed to be CP-even. But they contain terms 
proportional to sin0 g , which is CP-odd. This is accounted for because, in all cases, 
sin g is multipled by a term involving the helicity amplitudes which is also CP-odd. 
Similarly, cos0 g (CP-even) is multipled by a helicity-amplitude term which is also 
CP-even. The upshot is that the Ki(t) + Ki(t) (i = 1,2,3,5) are indeed CP-even. 
And it is straightforward to verify that the Ki(t) + K^t) (i = 4, 6) are CP-odd. 

The key point here is the following. The individual K^s and K^s [Eq. (j!4j) ] 
depend on four functions of time: e cos Amt, e sin Amt, e" ri cosh(Ar/2)t, 
and e _r * sinh (Ar/2)t. However, in the expressions above, the dependence on the 
functions e~ rt cos Amt and e~ rt sin Amt cancels, so that the untagged observables 
depend only on e~ rt cosh (AT/2)t and e~ rt sinh (AT/2)t. For B® mesons, Ar = 0, 
so that the untagged observables are equal to e _rt x simple sums of functions of 
the Ai and A\. On the other hand, since Ar ^ for B° s mesons, the untagged 
observables are now complicated functions of the Ai and A^. 

In addition, we have 



e~ rt cosh (AF/2)t = - (e- rii + e"^ 4 ) , e~ rt sinh (AT/2)t = ~ (e" rit - e"^*) . 



If the e _rit//2 and e~ THt ^ 2 terms can be distinguished experimentally, which is doable 
for decays, the untagged time-dependent angular distribution provides 12 observ- 
ables, 2 for each Ki(t) + Ki(t). Thus, B° s — > V1V2 decays are particularly interesting. 

2.4 Time-integrated untagged distribution 

As noted in the previous subsection, because Ar ^ for B® mesons, B° s decays can 
be treated without tagging. The time-integrated untagged angular distribution can 
be obtained by integrating the Ki(t) + Ki(t) observables over time: 



(17) 




(18) 



where 
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(Ki) = ^J dt{Ki(t) + Ki{t)) . 



(19) 



(Ki) 
(K2) 



One can obtain the (-fQ)'s from Eq. f)16p : 



2(1 -y 2 s ) 
2(1 -y 2 s ) 

2(1 -y s 2 ) 

2(1 -y 2 s ) 



\A \ 2 + \A \ 2 ) -2(Re(A*,4 o )cos0 s + Im(A*A ) sin0 s J?/, 
|A||| 2 + |A||| 2 ) -2(^Re(A|A||)cos0 s + Im(A|A||)sin^y, 
|-4±| 2 + |-4_l| 2 ) + 2(Re(A* ± A_ L ) cos0 s + Im^AjJ sin</> s Jy, 
Im(Aj.A:) - Im(^ ± A|)) - ((lm(A ± AD - Im(A|Aj!)) cos</ 



+ (Re(AiAJf) + Re(Aj_A|)) sin0 s W 



TB fl 



2(1 -yS 



Re(A|| Aq) + Re(A|| Aq) ) - (Re(A,|A*) + Re(A\\A* Q )) cos, 



<#6> 



2(1 



(Im(A||A*) -Im(A||A5))sin^)j/. 



Im(A ± A*) - Im(A ± A*; 



((Im(A x A*) - Im(Aj 



.^0))COS S 



+ (Re(A ± A*) + Re(A ± A*)) sin0 s W 



(20) 



where y s = Ar s /2r s . 

At this stage, one clearly sees the effect of a nonzero AT S (or y s ). For B® decays, 
ATd = 0, so there are no terms proportional to yd = ATd/2T ( i in the (Ki). Indeed, 
the (Ki) take the same form as the (Kit) + Ki(t))\ t =o [Eq. flBJ)]- However, this 
does not hold for 5° decays. Because of the nonzero y s , the (Ki), which are time- 
integrated quantities, take a different form than they did at t = 0. And this means 
that, if general b — > s NP is considered, the formulae relating certain observables to 
the (Ki) must necessarily include terms proportional to y s . As we will see, this holds 
specifically for the polarization fractions, CP-violating triple-product asymmetries, 
and the CP-conserving interference term. 
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2.5 Effective lifetime 



In general, the expressions for Ki(t) + Ki(t) [Eq. f fl6|) ] and (Ki) [Eq. (120]) ] have the 
form 



i^(t) + ATi(t) = 2e -lt A cosh (Ar/2)t + A sinh (AT/2)t 



r B s 



1 - 



A h + A V 



(21) 



where the experimental observables (dependent on the Ki) are on the left-hand side, 
and the theoretical expressions (dependent on Af 1 and A ) are on the right-hand 
side. (We have implicitly assumed that Ar ^ 0, which implies a 5° decay.) Af 1 and 
.A?' 1 can be related to the experimental observables via the effective lifetime [30] : 



T 



eff,i 



j Q ° C 't(Kj(t) + Kj(t))dt 



t Bs (1 + 2A i AT y s + y\ 



1 " 2/, 2 



;i+A^ 



where A\ r = Af/Af. 

Using Eqs. (EED and (J22), one can relate the ^4f to the (Ki): 



A 



- {Ki) (2 



eff,i 
B s 



^B s 



(i - y 2 s)) 



(22) 



(23) 



The A; h can be obtained from A 



Ar- 



3 B° s -> (jxj) - SM 

The results of the previous section are completely general. In this section we focus 
on the angular distribution of the pure b — > s penguin decay B° s — > 00 within the 
SM. 

In the SM, the amplitude for B® — > cfxp can be written 

a(b° s 00) = a«p; + \^P' C + \$P' U 

= + ^P'uc , (24) 

where A<f = V g ^V^ s . (As this is a 6 — > s transition, the diagrams are written with 
primes.) In the second line, we have used the unitarity of the Cabibbo-Kobayashi- 

Maskawa (CKM) matrix (A« +A^+A^ = 0) to eliminate the c-quark contribution: 
p/ — p/ p/ p/ — p/ p/ 

r tc — r t r ci r uc — r u r c- 
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Now, |A^| and |A« | are 0(A 2 ) and 0(A 4 ), respectively, where A = 0.23 is the 
sine of the Cabibbo angle. This suggests that the A^P„ C term can be neglected 
compared to A^ P/ c . However, if one does this, one must be consistent and neglect 
allO(X 4 ) terms. In particular, Im(A^ ) is 0(A 4 ), and so it too can be neglected. But 
since 2(3 S = — arg ((q/p)(A/ A)), one also has f3 s = because (q/p) = {A/ A) = 1 in 

(s) 

the limit where X t is real. Thus, in the approximation of neglecting all quantities 
of 0(A 4 ), there are no nonzero weak phases in B° s — > 00, either in the mixing or in 
the decay. 



3.1 Untagged distribution 

In the approximation of neglecting all weak phases in B® — > 
observables [Eq. f lTB]) ] are 



the untagged 



SM 



-17 



2|A)| 2 ( cosh(Ar/2)t - sinh (AT/2)t 



(K 2 (t)+K 2 (t)) S M = e 2|A||| 2 cosh (AT/2)t - sinh (AT/2)t 



(K 3 (t) + K 3 (t)) S M = e 



-rt 



2\A ± \ 2 [ cosh (AT/2)t + sinh (AT/2)t 



(^) + ^)W = o, 

(# 5 (*) + K 5 (t)) SM = e~ Tt 2Re(A,|A*) ( cosh (Ar/2)i - sinh (Ar/2)t 



(K 6 (t) + K 6 (t)) 



SM 







(25) 



We have Af 1 = =Fv4f fc [Eq. ( 12"T|) ]. where the minus sign is for z = 1, 2, 5, the plus sign 
for i = 3, and both quantities vanish when i — 4, 6. The effective lifetimes are then 
predicted to be 

r £s , ! r. r eff.i T B S 



'B S ,SM ~ i \ ' 5 5 5 'Bs,SM 



[i + y.) ' " ' (i .<//) ' ' (26) 

If the measurement of an effective lifetime differs from the SM prediction, this will 
be a sign for NP [30] . 

The SM untagged time-dependent angular distribution for B° s — > 00 takes the 
form 

d\v B ° + v B ') 



9 



.F L (g 2 ,u;)/C L (t) + r H {<f,w)K, H {t) 



dtdu 327T 
where the angular and time- dependent terms are 

F L {CS) = \\A \ 2 h{Q) + \A\\\ 2 f 2 {Q) + \A Q \\A\\\cos{5\\-5 )h{u) 



(27) 
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JC L (t) = 2e~ TLt = 2e~ r *(cosh(Ar/2)t-sinh(Ar/2)^ 

K H {t) = 2e~ rHt = 2e~ rt ( cosh (AT/2)t + sinh (AT/2)f) , (28) 

in which (S\\ - S ) = arg(A||Ao). 

Thus, if the e _rL */ 2 and e _rH */ 2 terms in the time-dependent angular distribution 
[see Eq. 017|) ] can be distinguished experimentally, the and cos — So) can be 
measured. However, as we will see in the next subsection, these observables can be 
obtained from time-integrated measurements. 

3.2 Untagged time- integrated distribution 

In the SM, the observables in the time-integrated untagged distribution are 

(K 1 ) = -^\A \ 2 , (K 2 ) = -^\A\\\ 2 



l + y s i + y, 

1-2/* 



(K s ) = \A X \ 2 , (K,) = 



(K s ) = -^-\A Q \\A\\\cos(6\\-6 Q ) , (tf 6 ) = 0. (29) 

We have y s = 0.088 ± 0.014 and r^ 1 = (0.6580 ± 0.0085) ps" 1 [30]. With this 
knowledge, the \ A^\ and cos (8n — So) can be extracted from the above measurements. 
This is what CDF and LHCb have presented 



3.3 Polarization Fractions 

With no weak phases in the decay, we have Ah = Ah, and the \Ah\ 2 can be measured 
in the untagged time-integrated distribution [Eq. ( 129]) ]. The polarization fractions 
are given by 







\A 


12 




\A 


l 2 + 


\ A \\ 


| 2 + 


l^l 2 






\ A \\ 


12 




\Ao 


l 2 + 


\ A \\ 


| 2 + 


l^l 2 






\A± 


2 




\A 


l 2 + 


(Ay 


P + 


\A±\ 2 



(30) 
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with total polarization f tot — fo + f\\ + f± — 1. 

Now, in the presence of NP the distribution changes, and so the experimental 
measurements have to be reinterpreted. We address this issue in the next section. 

4 B° s -> 00 - SM + NP 

In this section, we consider NP contributions to B° s — > <p<f), in the mixing and/or in 
the decay. 



4.1 Polarization Fractions 

The polarization fractions can be written as 

fo 

f\\ 
fx 









\A \ 2 + \Ao 


2 




A? 






2 + \ A 


2 + 


A \\\ 2 + \ A \\ 


2 + \A ± \ 2 + \A_ 


l| 2 


V A ch 


) 








\ A \\\ 2 + \ A \\ 
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Ach 
• Al 2 




\A 


2 + \ A 


2 + 


A \\\ 2 + \ A \\ 


2 + \A ± \ 2 + \A_ 


l| 2 


Z-a=l,2,3 


) 








A 1 _\ 2 + \A 1 


2 




Ach 


■ (31) 


\A 


2 + 1^0 


2 + 


A \\\ 2 + \ A \\ 


2 + \A ± \ 2 + \A_ 


l| 2 


V A ch 

Z-a=l,2,3 


above, 


the f h 


are 


written in 


terms of the 


A h \ 2 


and \A h \ 2 . 


However, 



as noted above, what is measured experimentally in the time-integrated untagged 
distribution are the {Kj). It is therefore necessary to express the fh in terms of the 
(Ki). This is done as follows. Using Eq. fl23|) . one can write 



•At = —0- + rhys)Yi, * = 1,2,3, (32) 
where the quantity Yi is related to T B or A l Ar : 

Yi = 7T-^ rf2- ^(1 - y 2 s )) = , (33) 



(l + Ws )V r Bs * S/ J (l + A\ T y s ) 
with 77i j2 = 1, and t] 3 = —1. From Eq. ( 120]) we have .4 4 A r = Af 1 j Af 1 = 

-2(Re(A*A Q )cos(j) s + lm(A*A )sm(p s ]/(^\A \ 2 + \A \ 2S ) , i = l, 
< -2fRe(A|A||)cos0 s + Im(A*A||)sin0 s j/nA||| 2 + |A||| 2 J , i = 2, (34) 
2(^Re(AlA ± )cos0 s + lm(A* L A ± )sm^ s J/(^\A ± \ 2 + \ A 



i = 3 
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In the SM, the weak phases of the Ah vanish and <fr s = 0, so that A\ r = ±1 (the 
minus sign is for % — 1, 2, and the plus sign is for % — 3). This implies that Y\,%z = 1, 
so that the polarization fractions are 



SM 



fo 



<^i)(i + y. 



(^i>(i+w.: 


) + (K 2 )(l + y s )- 




- 




(K 2 )(l + y s ) 








) + (K 2 )(l + y s )- 


f (^3>(l 


-y-) 












) + (K 2 )(l + y s )- 




-y s ) 



J± ' (K 1 )(l+y 8 ) + (K 2 )(l + y a ) + (K 3 )(l-y a ) ' 1 J 

Note that these are consistent with Eq. (1291) . However, if there is NP in the mixing 
and/or the decay, we have Yl,2,3 7^ 1> so that the polarization fractions take the form 





vy s . 


)Y 1 - 


-(K 2 )(l+y s 


)Y 2 H 














(K 2 ){l + y s . 


)Y 2 








<*1>(1- 


vy s . 


)Y 1 - 


-(K 2 )(l + y s [ 


)Y 2 H 




-y»! 










(K 3 )(l-y s . 


)Y 3 












)Y X - 


-(K 2 )(l + y s [ 


)Y 2 H 


-<Ka>(l 







,fo 

/ll 

/x = (Ki)(i + y s )yi + (K 2 }(i + L)r 2 + (K 3 )(i - y s )Y 3 • (36) 

The fh are expressed completely in terms of measured quantities. The (K^'s are 
obtained from the untagged angular distribution, and one can calculate the Y> using 
the measured effective lifetimes. If the effective lifetimes have not been measured 
then, v4Xr can be varied within a certain range to get a range for the Y$. 

Thus, to obtain the correct polarization fractions in the presence of NP, Eq. fl36|) . 
which includes factors of Yj, must be used. This is one of the main points of the 
paper. However, experiments have used Eq. fl35|) . so they have effectively excluded 
NP. If this possibility is allowed, the analysis must be redone and we discuss this in 
Sec. 5. 

The difference between Eqs. (I3"5|) and (I3T)|) is related to the difference Y^ — 1. One 
can see from Eq. ( 13"B"|) that Y^ — 1 — > in the limit that y s — > 0. This indicates that 
fh — fh M = 0(y s ). Since y s = 0.088 ± 0.014, this corresponds to a correction to the 
polarization fractions of O(10%). This is not large, but it may be important given 
that the measurements hope to identify the presence of NP. 



4.2 Other Observables 

In Sec. 12.2} we noted that the angular distribution of the decay B® — > V\V 2 (q = 
d,s) is proportional to J2i=i Ki{t)fi{oj), where u = (cos 9%, cos 9 2 , $) [Eq. (jBJ)]. In 
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the previous subsection, we discussed polarization fractions, observables which are 

dependent on (Ki), i — 1, 2, 3. We now turn to i — 4, 6. 

In the present case, K± and K$ are related to the triple products (TP's) in 
— > cjxj). The expressions for the untagged observables in B® s — > V1V2 are given 

in Eq. f lT6|) . For convenience, Ki(t) + Ki(t) (i = 4,6) are repeated below: 

K A (t) + K A (t) = e- Fat \(lm(A ± AD - Im(A ± A|)) cosh (Ar s /2)t 



((Im(A ± Af) - Im(A ± Af)) cos. 



+ (Re(Aj.A:) + Re(Aj.A:)) sin0J sinh (AT s /2)t 



K e (t) + K 6 (t) = e~ Vst [(lm(A ± A*) - lm(A ± A* )^ cosh (AF s /2)t 
- ((Im(A ± A*) - lm(A ± A* )) cos0 s 



+ (Re(A ± A*) + Re(A ± A* )) sin <p s ) sinh (AT s /2)t 



(37) 



Now, as discussed earlier, in the SM the weak phases in — > (fxp, both in the mixing 
and in the decay, are all approximately zero, so that K±(t) + K^{t) and K$(f) -\-K${£) 
vanish. Thus, if one finds evidence for a nonzero TP, this is a clear sign of NP. 

Suppose first that there is NP, with a nonzero weak phase, only in the mixing. In 
this case, the first two terms in each of Ki(t) + Ki(t) (i = 4, 6) are zero, but the third 
is nonzero. This is a particularly interesting situation, as it corresponds to a TP 
generated through mixing. It arises only because Ar s is nonzero; mixing-induced 
TP's cannot be produced in B® decays. And, although Ar s 7^ 0, it is still not large, 
so that the associated TP is also rather small. 

The second possibility is that there is NP, with a nonzero weak phase, only in 
the decay. In this case, the first two terms in each of Ki(t) + Ki(t) (i = 4,6), 
proportional to cosh(Ar s /2)t and cos0 s = 1, are nonzero, but the third is zero. 
And of course one can have NP in both the mixing and the decay. If a TP is seen, 
its source can be determined through its time dependence. 

Both Ki(t) + Ki(t) {i = 4, 6) are CP- violating, so they correspond to true TP's. 
They can be nonzero only if there are two interfering amplitudes with a relative 



weak phase. If there is NP in the mixing, the amplitudes are A(B S 



A(B° -+ Bl 



; if there is NP in the decay, the amplitudes are A(B® — > 



and 



SM 



and A{B° S — y 4>4>)np- In addition, in order to produce a TP, the two interfering 
amplitudes must be kinematically different [24] . For the case of NP in the decay, 
this is satisfied straightforwardly. But for NP in the mixing, how are B® — > cfxp and 
B® — y B® — y cfxp kinematically different? The point is that mixing-induced TP's are 
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generated due to a nonzero AT S . That is, although B® 00 is a penguin decay, 
B° s — > <jxj) occurs via a mechanism which contributes to Ar s . For example, one 
possibility is the B° s — > B° s transition via the intermediate states D* + D*~ [34] , with 
the Bg decaying to (fxfi. The B° s and B Q S decays are clearly kinematically different. 

We now turn to the measurement of TP's. Here we focus on the time-integrated 
untagged observables, (Ki). We have (Ki) oc Af 1 + Af h y s [Eq. fl2Tj) ]. Specifically, 
the (A 4j6 ) are given in Eq. (120]) : 



(K 4 ) 



2(1 - y* 



hn(A ± Al) - lm{A ± A: 



(lm(A ± AT l ) - lm(A ± At)) cos. 



T B S 



2(1 - y\ 



+ (Re(A ± A|) + Re(A ± 4))sin0 s? /, 



(lm(A x A*) - lm(A ± A*)^j - ((Im(A ± A*) - lm(A ± A*)) cos, 



+ (Re{A ± A* ) +Re{A ± A*))sm<f> s )y 



(31 



The TP's in the untagged distribution can be measured by constructing asymmetries 
involving the angular variables. We start by integrating Eq. ( TI~8l) over 9\ and 02 to 
obtain the differential rate: 



d(T(B° q -» V 1 V 2 )) _ ^ 
d$ 2?r L 



(i^i) + 2(A 2 ) cos 2 $ + 2(A 3 ) sin 2 $ - 2(A 4 ) sin 2$ (.39) 



Note that the time-integrated untagged decay rate can be obtained by integrating 
out the azimuthal angle 



(r(R° -> v^)) = (^i) + (a 2 > + (a 3 ) 



(40) 



Following Ref. j2U] we can define asymmetries to measure the TP's. We begin 
with i = 4, for which f^{uj) = —2 sin 2 61 sin 2 # 2 sin 2$. We define u = sin 2$. The 
TP asymmetry between the number of decays involving positive and negative values 
of u is given by [201 124] 



1 v (T(B° 8 -)■ #), » > 0) - gxgj ^ < 0) 

u 2 L <r(s° -> #), u > 0) + (r(s° o) 



-L4 (2) l 



TT 



U (2) l 
ly^r J ex p 



(*4> 



(41) 



As noted above, if A u ^ is found, this would clearly indicate NP. However, we 
would like to know the relation between [^l^jexp and the theoretical expression for 
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the TP in Eq. (138|) . The measured TP [A^ \ ex p is related to [A? ]theo y i a 



r.(2), _ r .(2), t_b^ 

\exp IrVT \theo ^p^O ^ 



where A { £ T = Af / Af and 



[A 



(2)1 

7 1 Jt/ieo 



■A 



1 



(lm.{A L AV) - Im(AxAn)) 



(42) 



(43) 



If we define the dimensionless theoretical TP as 

TV 2 = [Ar'jtheo 



(2) ] TB„ 



(44) 



Eq. ( 142]) details the corrections to the naive relation [A^ ] eX p — TV 2 due to a nonzero 
(NP) Af. (In the SM, [AP] theo = 0, so the relation is trivial.) 

For i = 6we have fe(uj) = — v / 2sin26' 1 sin 202 sin $. We define 
v = sign(cos0i cos #2) sin$, which has the following associated TP asymmetry [20J: 

A .irm ^ 



v > 0) - (T(B° S 00), v < 0) 

(r(s° 00),^ > 0) + (r(s° 00), u < 0} 



Then 



where A^ r = Af / Af and 



-— U (1) l 



7T 



[A 



(Dl 

T 1 exp 



[A. 



(1)1 _ r 4(1)1 T Bs 



(r(5 s ° 



1 + A&y. 



(45) 



(46) 



[Asp ^]theo — A§ 



(lm(A ± A* ) - lm(A ± A*)) 



(47) 



We can again define the dimensionless theoretical TP as 



TVl = [Asp]theo 



(r(5 a o 



Eq. ( 146]) gives the corrections to the naive relation [A T 



(1)1 



J cxp 



TPi. 



(4* 



Finally, we turn to z = 5, which corresponds to a CP-conserving observable. 
From Eq. flU}, 



2(l-2/ s 2 



Re(A||A*) + Re(A||A*; 
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(Re(A||A*) + Re(A||A*)) cos< 



(Im(A||A*) - Im(A|,A*)) sm<j) s )y, 



(49) 



We have f$(uj) = y/2 sin 26i sin 29 2 cos $, so we define u> = sign(cos9i cos9 2 ) cos $ 
The associated asymmetry is 



We have 



A 



gXg -» > 0) - (T(gg -» w < 0) 

(T(B° -)> 00), w > 0) + (r(5° #),«;< 0) 

[A( 5 )] e:cp =(K 5 )/(r(5 s ^ 



7T 



;i + A®y s 



r e//,5 



(50) 



where = A^/A^ 1 , the effective lifetime t^' 5 is defined in Eq. ( |22l) . and 
[A (5) Wo = = ^(Re(A {l A* ) + Re(A||AS; 



(51) 



(52) 



4.3 NP Parameters 

12 observables can be measured from the time-dependent untagged angular distri- 
bution (Sec. 12. 31) . With these, one can identify if NP is present in the mixing and/or 
the decay. However, we will also want to identify its properties. To be specific, if 
there is NP in the decay amplitude, it will be important to measure the various NP 
parameters. With this in mind, the question is: how many theoretical unknowns are 
there in the most general SM + NP B° s — > (fxfi amplitude? If there are fewer than 
12, then we can extract all the unknowns. 

In writing the SM + NP — > (jxj) amplitude, we have the following points: 

• The SM weak phases are ~ 0. 

• Assuming that the NP amplitudes satisfy |*A^ P | < |.4.f M |, the NP strong 
phases are negligible [31]. This means that if there are many NP amplitudes, 
they can all be combined into a single term with an effective magnitude and 
weak phase. 

• In the heavy-quark limit, we have A^_ M = — Afi M [32] . 
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Taking these points into account, the most general SM + NP B® — > <jxp helicity 
amplitude can then be written 

Ah = \A s h M \e< M + K p |e #h • (53) 

There are a total of 11 unknown theoretical parameters - 5 magnitudes (2 SM, 3 
NP), 2 SM strong phases, 3 NP weak phases, and the mixing phase <p s . In principle, 
these can all be extracted from the 12 observables. 

However, note that Eq. ( 153]) includes a different NP weak phase <ph for each 
helicity amplitude. But in many NP models the weak phases are helicity indepen- 
dent. In this case there is only one NP weak phase 0, and the number of theoretical 
unknowns is reduced to 9. This is a model-dependent result, but it is still very 
general. 

Finally, if the NP is purely left-handed or right-handed, then = T^j| P |33j . 
which further reduces the number of theoretical unknowns by one. 

In all cases, assuming the time-dependent untagged angular distribution can be 
measured, there are more observables than unknowns, and so we will be able to 
extract all the NP parameters in the decay. In this way, we may be able to identify 
the type of NP that is present. 

5 Numerical Analysis 

Recently, the CDF and LHCb Collaborations have reported measurements for the 
polarization amplitudes, the strong-phase difference between A\\ and Aq, and the 
triple-product asymmetries in B° s — > cjxf). The LHCb results [22] are summarized 
in Table [TJ The values are in good agreement with those reported by the CDF 
Collaboration [2TJ, except for the TP's, though all measurements are consistent 
within errors. 



Observable 


Measurement 


1 A I 2 


0.365 ± 0.022 (stat) ± 0.012 (syst) 


1 A . |2 
1 -L 1 exp 


0.291 ± 0.024 (stat) ± 0.010 (syst) 


\A..\2 
1 II \exp 


0.344 ± 0.024 (stat) ± 0.014 (syst) 


cos(<5|| - S ) 


-0.844 ± 0.068 (stat) ± 0.029 (syst) 




-0.055 ± 0.036 (stat) ± 0.018 (syst) 


A 


0.010 ± 0.036 (stat) ± 0.018 (syst) 



Table 1: Measured polarization amplitudes, strong-phase difference, and triple- 
product asymmetries in B® — > <p(f) [22J. The sum of the \Ah\ 2 exp terms is constrained 
to unity. 

The experiments have measured the (Ki) and constructed the polarization frac- 
tions assuming the SM. As discussed previously, if one allows for the possibility of 
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NP in b — y s transitions, this analysis must be modified. This is done below. 
We denote the measured value of y s as y S Q. From Eq. (1291) we have 





\M 2 = 




i + y s 


1 + VsO 




\M 2 = 


TB S 


i + y s 


1 + VsO 


r B s 


\A ± \ 2 = 


TB s 


1 - Vs 


1 - VsO 



(Ki) = T ^\A \ 2 = -^\A 

l + Vs 1 + VsO 

(K2) = ^\A ll \ 2 = -^\A ll 



2 

ys=y s o ' 



2 

ys=y B o ' 



l^lUo • (54) 
The experimental measurements in Table [1] are then 



\A 



|2 _ I » | y a =y s p 



|2 1 IMj/3=3/sO 



\A±\ 2 - 

\A±\ exp = , A , 9 . i ^ i9 ~T1 i9 • (55) 



Vs=y s o 



One can now calculate the polarization fractions in the SM as a function of y s . 
Inputting the expressions for the (Ki) from Eq. ( 1541) into Eq. ( 135|) . and using Eq. ( 1551) . 
we obtain 



Jo 



14)1 



2 i I /I ,.12 i I 4 12 



l^-olLpi+^ + l^nlLpi+^ + 1^ 

i+y. 



fi 



SM _ 1 H |ea; Pl+?/ s o 



\A\\\ 2 exf 



MJ2 1 +^ I 1/1 1, 12 l±Vs_ , I 4 12 



I A 12 

/f M = tt 1 ±lexpl ~y°° — . (56) 



Hence the |A|L P in Tabled] are just the ff M defined in Eq. (1561 with y s = y s0 . 

The true polarization fractions can then be obtained by inputting the expressions 
for the (Ki) from Eq. (ISlfl) into Eq. (|36|) . and using Eq. (|55|) : 



/o 



I A 12 

I^O| e xpT+^-^l 
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/II 



12 1+Vs 
1 e-xpl+y s0 



4J2 1+ y s V, i I A,, 12 V-, i 1 ,4 

"^Oleaipl+jfco 1 1 + I^IILttiij-,. „ r 2 "t 



12 i-ya 



/J 



1- 



Y, 



AJ 2 1+Va Y, 4- I Au I 2 



ea 'P 1+VsO 



Y 2 + \A_ 



12 i-ys 
\expl-y s0 



Y. 



(57) 



In Fig. [T] we plot the dependence of the polarization fractions / , /y and /j_ as a 
function of y s . This figure is read as follows. In all plots the horizontal region 



represents the experimental result, in which 



j ex p 



is allowed to vary by 



±lcr (see Table [T]). Also, the vertical bands correspond to y s , with ±lcr (green) 
or ±3a (yellow) errors. In the SM we have Yi = 1, corresponding to {A Ar = 
-1,A 2 AT = -1, A 3 AT = 1) [Eq. ([M])]. In order to illustrate the effect of NP, we take 
(A l Ar = l,A 2 Ar = -l,A 3 Ar = -1) (red line) or (AX r = -1, A 2 AT = 1,A 2 AT = 1) (blue 
line). For these values of A Ar , we have Yi ^ 1. Consider first Jq. In the SM the 
experimental measurement implies 0.33 < fo < 0.40. However, with NP, the value 
of fo can lie outside this range - for example, on the red line it can be as small 
as 0.29. The behavior is similar for fa and f±. This shows explicitly that, in the 
presence of NP, the B° s — > <jxf) polarization fractions can be changed from their SM 
values by 0(10%) for the current value of y s . 
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Figure 1: The dependence of the theoretical polarization fractions fo, fa and f± on 



the decay width parameter y s for different values of -A^p 3 - The red line corresponds 



to (A Ar = l,A Ar = —l,A Ar = — 1), while the blue line has (A AT 



■ 1 A 2 
J-; 71 Ar 



1 Ar — J-j^Ar — J-i^Ar - 
l,A 2 Ar = 1). In all plots the experimental result [|A|^ =0 u ± ] exp (horizontal region) 

is allowed to vary by ±ler (see Tabled]). The vertical bands correspond to y s , with 
±1(7 (green) or ±3cr (yellow) errors. 



The relation between A u and [A}p ]theo is given in Eqs. ( )4T|) and ( 142]) ; that be- 
tween Av and [A^ \theo is given in Eqs. ( 145]) and (146]) . These can be rewritten 



as 



r,(2)i r Ba 

i T \theo rp(B° -»■ 



2 d + 4 4 k 
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Figure 2: The dependence of the theoretical TP's [AP} theo T B J (T(B° S -> 00)) (left) 
and [At ]theoT~B s / (r(-B° — > 00)) (right) on y s for different values of A^p . In the red 
and blue regions, we take A^r = ±lj respectively. Also, .Au (left) and A v (right) 
are allowed to vary by ±ler (see Table [T]). The dashed black lines correspond to 
the central values of A u (left) and A v (right) with A^jj) = 0. The vertical bands 
correspond to y s , with ±lcx (green) or ±3<r (yellow) errors. 



,2"> 



r .(i) 1 r Bs _ 7T u-yj , . 

^ ^ (r(fi? -> 00)) - y/2 (1 + A Wy a) ■ (58) 

In Fig. [2] we plot the dependence of the theoretical TP's [A^]^ eo tbJ (r(B® — >■ 00)) 
and [A^]^ eo TB s /(r(S° — > 00)) as a function of y s . The dashed black lines corre- 
spond to the central values of A u (left) and A v (right) with A^j} = 0. The vertical 
bands correspond to y s , with ±lcr (green) or ±3cr (yellow) errors. In the red and 
blue regions, we take A^ v = ±1, respectively, and allow A u (left) and A v (right) 
to vary by ±1ct (see Table [T|). It is clear from these figures that, in the presence 
of NP, the values of the theoretical TP's can differ significantly from the measured 
asymmetries. (This is not surprising since the TP's vanish in the SM.) 

Finally, for i = 5, we have estimated the measured value of the CP-conserving 
observable as follows: 

[A ( % xp = \A \ exp | AhI^ cos(5|| - 6 ) = -0.299 ± 0.030 . (59) 

The relation between [v4^ 5 ^] e:r p and [A( 5 )] t /j eo is given by [see Eq. fl5"Tj) ] 



[A{%h - (rwS 00)) " [Ai%xp (i + aIL) • (60) 

In Fig. [3] we plot the dependence of [A^] t heoTB s / (r(J3° — > 00)) as a function of y s . 
As before, the value of this quantity can differ from Eq. fl59l) by as much as O(10%) 
for the current value of y s . 
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Figure 3: The dependence of [A ( T '} theo T B J (T(B® — > 00)) on y s for different values 
of A\ v . In the red and blue regions, we take A\ T = ±1, respectively. Also, [A^ 5 ^]^ 
is allowed to vary by ±lcr [see Eq. (|59|) ]. The dashed black line corresponds to the 
central value of [A^]^ with A 5 Ar = 0. The vertical bands correspond to y s , with 
±l(j (green) or ±3<r (yellow) errors. 



6 Conclusions 

The main goal of studying the B system is to find evidence for physics beyond the 
standard model (SM). One possibility is new physics (NP) in b — > s transitions. At 
present its status is uncertain. It seems unlikely that the effect of such NP can be 
very large, but a smaller effect is still possible. In this paper, we consider b — y s 
NP. However, in contrast to what is usually done, i.e. considering only NP in B®-B® 
mixing, here we also allow NP in the decay. In particular, we examine the effect of 
such NP on the angular distribution of B® — > V1V2 (<7 = d, s), where V\ t 2 are vector 
mesons. 

Our principal result is the following. The parameters of the untagged, time- 
integrated angular distribution can be measured experimentally, and certain ob- 
servables can be derived from these parameters. However, in the presence of NP, 
the formulae which relate the parameters to the observables must be modified from 
their SM forms. We find six observables for which the relation between the experi- 
mental data and theoretical parameters must be modified, corresponding to the six 
terms (i = 1-6) in the angular distribution. For i = 1-3 they are the polarization 
fractions, for i = 4,6 they are the CP- violating triple-product asymmetries, and i = 5 
corresponds to a CP-conserving observable. The modifications for the polarization 
fractions are most interesting. These are due in part to the nonzero width difference 
in the B®-B® system, and so are important only for B® decays. In particular, there 
can be important effects on the pure b — > s penguin decay B° s — > 00. 

In light of this, we re-analyze the B° a — > 00 data to see the effect of these 
modifications. AT S /2T S ~ 10%, so that the modifications of the formulae lead to 
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0(10%) changes in the polarization fractions. These are not large, but may be 
important given that one is looking for signals of NP. 

Finally, if the NP contributes to the b — > s decay, we show that the measurement 
of the untagged time- dependent angular distribution provides enough information - 
12 observables - to extract all the NP parameters. 

Note added: while the paper was being written, D0 produced a direct measure- 
ment of the semileptonic charge asymmetry in B° s decays [36] , and they say that 
it agrees with the SM. Technically, this is true. The SM predicts a s sl ~ 2 x 10 -5 , 
and D0 measures a s sl = [-1.08 ± 0.72 (stat) ± 0.17 (syst)]%. While this result is 
consistent with 0, the errors are large enough that NP is also a possibility. 
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